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We study flux tubes on Higgs branches with curved geometry in super- 
symmetric gauge theories. As a first example we consider 7V= 1 QED with 
Qh! one flavor of charges and with Higgs branch curved by adding a Fayet- 

Iliopoulos (FI) term. We show that in a generic vacuum on the Higgs 
branch flux tubes exist but become "thick". Their internal structure in 



the plane orthogonal to the string is determined by "BPS core" formed by 
heavy fields and long range "tail " associated with light fields living on the 
Higgs branch. The string tension is given by the tension of "BPS core" 
plus contribution coming from the "tail". Next we consider A/ = 2 QCD with 
gauge group SU{2) and Nf = 2 flavors of fundamental matter (quarks) 
with the same mass. We perturb this theory by the mass term for the 
adjoint field which to the leading order in perturbation parameter do not 
break M= 2 supersymmetry and reduces to FI term. The Higgs branch has 
Eguchi-Hanson geometry. We work out string solution in the generic vac- 
uum on the Higgs branch and calculate its string tension. We also discuss 
if these strings can turn into semilocal strings, the possibility related to the 
confinement/deconfinement phase transition. 



1 Introduction 



The mechanism of confinement as a dual Meissner effect arising upon conden- 
sation of monopoles was suggested a long ago [JJ. Once monopoles condense 
the electric flux is confined in the (dual) Abrikosov-Nielsen-Olesen (ANO) 
flux tube [21 IS] connecting heavy trial charge and anti-charge. The flux tube 
has constant energy per unit length (string tension ). This ensures that the 
confining potential between heavy charge and anti-charge increases linearly 
with their separation. However, because the dynamics of monopoles is hard 
to control in strong coupling gauge theories no quantitative description of 
this phenomenon was constructed in QCD. 

The revival of interest to this problem occurs after the work of Seiberg and 
Witten jUE]. They showed that in supersymmetric gauge theories the holo- 
morphy and electromagnetic duality are powerful tools to study non-Abelian 
dynamics at strong coupling. Using the exact solution of N = 2 supersym- 
metric gauge theory they were able to demonstrate that the condensation of 
monopoles do really occurs near the monopole vacuum once N = 2 theory is 
slightly perturbed by the mass term of the adjoint matter [3j. 

Since then a lot of papers study confinement and formation of flux tubes 
in supersymmetric gauge theories |H1 IZl IH1 13 UHl HH [T21 H31 [TBI - 

One important particular feature of supersymmetric theories is the presence 
of modular spaces - manifolds on which scalar fields can develop arbitrary 
VEV's. If the gauge group is broken down completely by scalar VEV's such 
vacuum manifolds are called Higgs branches. In this paper we study the for- 
mation of flux tubes in supersymmetric gauge theories with Higgs branches. 

In fact , this problem was addressed earlier in ref . [TO] for the case of Higgs 
branches with flat geometry. In particularly, the formation of flux tubes 
on the Higgs branch of Af= 2 QCD with gauge group SU (2) and Nf — 2 
flavors of fundamental matter hypermultiplets (quarks) with the same mass 
was studied. The Higgs branch in this case represent an extreme type I 
superconductor with vanishing Higgs mass. It was shown in JU] that the 
flux tubes still exist in this set up although becomes logarithmically "thick" 
due to the presence of light scalar field. Because of this the confining potential 
behaves as 

with the separation L between two heavy well-separated (magnetic) charges. 
This potential is still confining but is no longer linear in L. Note, that in this 
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case we have condensation of scalar components of quarks (electric charges) so 
these are monopoles which are confined. As the potential between heavy trial 
charges is an order parameter to distinguish between different phases of the 
theory, we see that we have a new confining phase with a non-linear potential 
which is specific for Higgs branches in supersymmetric gauge theories. 

In this paper we continue to study flux tubes on Higgs branches concen- 
trating on the case of Higgs branches with curved geometry. The curvature 
on a Higgs branch is induced by adding of the Fayet-Iliopoulos (FI) term to 
the theory. 

Our first example is M= 1 QED with one flavor of quarks and the FI D- 
term added. This theory has a two dimensional Higgs branch of a hyperbolic 
form. In order to regularize this theory in the infrared one would like to lift 
the Higgs branch making scalars which fluctuate along the vacuum manifold 
slightly massive. It turns out that this can be done without destroying of 
M— 1 supersymmetry. We start with Af= 2 QED with one quark flavor which 
has no Higgs branch. Then we add a mass term for the neutral matter field 
breaking A/ = 2 supersymmetry down to M= 1 . Then in the limit of mass 
parameter of this field /i->oowe recover M= 1 QED with its Higgs branch. 
At large but finite \i we have Higgs branch lifted but the potential along the 
would be Higgs branch becoming more and more flat as we increase //. 

In this set up we study the formation of flux tubes in a generic vacuum 
on the would be Higgs branch. We show that flux tubes exist and have 
the following structure in the plane orthogonal to the string axis. They are 
formed by "BPS core" which size is determined by the Fl-parameter and 
long range "tail" formed by light scalar fields. The string tension is given by 
the tension of "BPS core" plus contribution coming from the "tail". If the 
VEV of scalar fields squared is of the same order as the Fl-parameter the 
contribution of the "tail" becomes a small correction to the tension coming 
from the "BPS core". That is why we can call this string "almost BPS". 
Note however, that this string is not a BPS state and belongs to a long 
supermultiplet. 

This string gives rise to the confinement of monopoles with the potential 
which depends on the monopole separation L as 




(1.2) 



Here the first term comes from the "BPS core" while the second one is 
determined by the "tail" . 
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Next we consider M=2 QCD with gauge group SU(2) and two flavors of 
quarks with the same mass. This theory has four dimensional Higgs branch 
jS]. We perturb this theory by the mass term for the adjoint field which to the 
leading order in perturbation parameter do not break M= 2 supersymmetry 
and reduces to FI term [TJE]- The non-zero Fl-parameter induces curvature 
on the Higgs branch. It becomes a hyper-Kahler manifold with Eguchi- 
Hanson ^7j geometry. 

We work out the ANO string solution in a generic point on the Higgs 
branch and calculate the string tension. Similar to the M= 1 case it is given 
by the contribution of "BPS core" plus correction coming from the "tail" of 
light fields living on the Higgs branch. 

We also discuss if this ANO string can turn into semilocal string ( see [TH] 
for a review on semilocal strings) by increasing the size of its "BPS core". 
If this happen this would ruin confinement making the potential between 
monopoles to fall-off as a power of L at large L. We show that if VEV of 
the quark fields is much larger then the Fl-parameter the semilocal string is 
not developed. ANO string appears to be stable so monopoles are confined 
by the potential (jl.2j) . 

However, for vacua at the base S2 cycle of the Eguchi-Hanson manifold 
for which quark VEV is equal to the Fl-parameter ANO strings turn into 
semilocal strings P^ ll6j. For these vacua we have deconfinement phase. We 
put forward a conjecture that confinement /deconfinement phase transition 
occurs exactly at the base cycle. Finally we consider effects which break 
M— 2 supersymmetry down to M= 1 and show that in this case the decon- 
finement phase disappears. 

The paper is organized as follows. In Sect. 2 we review flux tubes on 
Higgs branches with flat geometry [TO] . In Sect. 3 we introduce Af— 1 QED 
with FI term and in Sect. 4 study flux tubes on the Higgs branch in this 
theory. In Sect. 5 we consider Af=2 QCD with gauge group SU(2) and 
two flavors of quarks with the same mass and study flux tubes on the Higgs 
branch in this theory. In Sect. 6 we discuss the issue of semilocal strings. 
Sect. 7 contains our brief conclusions. 

2 Extreme type I strings 

In this section we review classical solution for ANO vortices in theories with 
flat Higgs potential which arises in supersymmetric settings [TU] . In particu- 
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lar, in ref. ^U] flux tubes on the Higgs branch of Af=2 QCD with two flavors 
of fundamental matter (quarks) were studied. Consider the Abelian Higgs 
model 

Sah = Jd"x^F^ + \V,q\ 2 + X(\q\ 2 -v 2 rY (2.3) 

for the single complex field q with quartic coupling A = 0. Here V M = 
<9 M — in e An, where n e is the electric charge of the field q. Following we 
consider first this model with small A, Then we take the limit A —>■ 0. 

The field q develops VEV, q = v , breaking down the U(l) gauge group. 
Photon acquires the mass 

m 2 = 2n 2 e g 2 v 2 , (2.4) 
while the Higgs mass is equal to 

m\ = AXv 2 . (2.5) 

The model ()2.3|) is the standard Abelian-Higgs model which admits ANO 
strings [3 Ej- For an arbitrary A the Higgs mass differs from that of the 
photon. The ratio of the photon mass to the Higgs mass is an important 
parameter, in the theory of superconductivity it characterizes the type of 
superconductor. Namely, for m q < m 1 we have the type I superconductor 
in which two well separated ANO strings attract each other. Instead for 
m q > rrij we have the type II one in which two well separated strings repel 
each other. This is related to the fact that scalar field produces an attraction 
for two vortices, while the electromagnetic field produces a repulsion. 

Now we consider the extreme type I limit in which 

m q <C m 7 . (2.6) 

We also assume the week coupling regime in the model ()2.3|) A <C g 2 <C 1 . 

The general idea to find the string solution is to separate behavior of 
different fields at different scales present in the problem due to the condition 
(12. This method goes back to the paper by Abrikosov [2] in which the 
tension of type II string under condition m q ^> m 7 has been calculated. The 
similar idea was used in JU] to calculate the tension of the type I string under 
condition m q -C m 7 . 

To the leading order in log m 7 /m g the vortex solution has the following 
structure in the plane orthogonal to the string axis ^U] • The electromagnetic 
field is confined in a core with the radius 

1 , m~ 

R log-^ . (2.7) 

m 7 m q 
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The scalar field is close to zero inside the core. Instead, outside the core, the 
electromagnetic field is vanishingly small while the scalar field behavies as 



\ logl /m q R g ) 

where r and a are polar coordinates in the plane orthogonal to he string 
axis. Here K is the Bessel function with exponential fall-off at infinity and 
logarithmic behavior at small arguments, Kq(x) ~ logl/x at x — > 0. The 
reason for this behavior is that in the absence of the electromagnetic field 
outside the core the scalar field satisfies a free equation of motion and (|2.8|) 
is a solution to this equation. From ()2.8|) we see that the scalar field slowly 
(logarithmically) approaches its boundary value v. 
The result for the string tension is ^U] 

2nv 2 , . 

T = — . 2.9 

log m^/rriq 

The main contribution to the tension in (J2.9|) comes from the logarithmic 
"tail" of the scalar q. It is given by the kinetic term for the scalar field in 
(12. This term contains a logarithmic integral over r. Other terms in the 
action are suppressed by powers of log m 7 /m g as compared with the one in 

(USD- 

The results in f!2.7j) . ()2.9j) mean that if we naively take the limit m q — > 
the string becomes infinitely thick and its tension goes to zero [TOj- This 
means that there are no strings in the limit m q = 0. The absence of ANO 
strings in theories with flat Higgs potential was first noticed in [T9*] . 

One might think that the absence of ANO strings means that there is 
no confinement of monopoles in theories with Higgs branches. As we will 
see now this is not the case |T0] . So far we have considered infinitely long 
ANO strings. However the setup for the confinement problem is slightly 
different. We have to consider monopole-anti-monopole pair at large but 
finite separation L. Our aim is to take the limit m q — > 0. To do so let us 
consider ANO string of the finite length L within the region 

— < L < — . (2.10) 

m 7 m q 

Then it turns out that 1/L plays the role of the 7i?-cutoff in Eqs. ()2.7p and 
(12. 9j) instead of m q jTOj- The reason for this is that for r ^ L the problem 
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is two dimensional and the solution of the two dimensional free equation of 
motion for a scalar field is given by (|2.8|) . If we naively put m q = in this 
solution the Bessel function reduces to the logarithmic function which cannot 
reach a finite boundary value at infinity. Thus as we mentioned above the 
infinitely long flux tubes do not exist. This was noticed in [T^j. However for 
r ^> L the problem becomes three dimensional. The solution for the three 
dimensional free scalar equation of motion behaves as g— v ~ l/M; where x n , 
n = 1, 2, 3 is a coordinate in the three dimensional space. We see that with 
this behavior the scalar field reaches its boundary value at infinity. Clearly 
1/L plays a role of IR cutoff for the logarithmic behavior of the scalar field. 

Now we can safely put m q = 0. The result for the electromagnetic core 
of the vortex becomes 

R g — logm 7 L , (2.11) 



m 7 



while its string tension is given by [TU] 

We see that the ANO string becomes "thick" but still its transverse size 
R g is much less than its length L, R g <L. As a result the potential between 
heavy well separated monopole and anti-monopole is still confining but is no 
longer linear in L. It behaves as ^Uj 

V(L) = 2nv 2 - — — - . (2.13) 
log m 7 L 

As we already explained in the Introduction the potential V(L) is an order 
parameter which distinguishes different phases of a theory (see, for example, 
review [20 ). We conclude that we have a new confining phase on the Higgs 
branches. It is clear that this phase can arise only in supersymmetric theories 
because we do not have Higgs branches without supersymmetry. 

3 J\f= 1 QED with Fayet-Iliopoulos term 
3.1 The model 

The field content of M= 1 QED is as follows. The vector multiplet contains 
U(l) gauge field and Weyl fermion A Q , a = 1,2. The chiral matter 
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multiplet contains two complex scalar fields q and q as well as two complex 
Weyl fermions ip a and ip a . The bosonic part of the action reads 

Sqed = J d A x 1^2*2, + V M gV M g + V M gV„g + g)| , (3.1) 

where = — V M = <9 M + f-A^, so we assume that that matter fields 
have electric charge n e = 1/2. The potential of this theory comes from the 
.D-term and given by 

V(q,q) = 9 j(\q\ 2 -\q\ 2 -^)\ (3.2) 

where parameter £3 arises if we include FI .D-term in our theory. We denote 
here the FI .D-term parameter £3. This notation will become clear later in 
7V= 2 setup where D and F FI parameters form a triplet with respect to 
global SU{2) R . 

The vacuum manifold of the theory (|3.1|) is a Higgs branch determined 
by the condition 

|<g>| 2 -|<g>| 2 = £ 3 . (3.3) 

The dimension of this Higgs branch is two. To see this note, that we have 
two complex scalars (four real variables) subject to one constraint (|3.3|) . Also 
we have to subtract one gauge phase, thus we have 4-1-1=2. 

Now let us consider the mass spectrum of the QED (|3.1|) . As it is clear 
from (j3.3J) at any non-zero £3 scalar fields develop VEV's breaking U(l) gauge 
symmetry. The photon mass is given by 

2 ^ 2 2 
= -g v , 

where we introduce the VEV of scalar field 

v 2 = |< q >| 2 + |< q >| 2 . (3.5) 

To find matter masses we diagonalize the 4x4 mass matrix for scalar 
fields in (13. 2 j) . It has three zero eigenvalues one of which is "eaten" by the 
Higgs mechanism while two others correspond to chiral massless multiplet 
living on the Higgs branch. The remaining fourth eigenvalue is equal to the 
photon mass ()3.4D . 

m H = m 7 . (3.6) 
7 



(3.4) 



This eigenvalue corresponds to the scalar superpartner of the photon in the 
massive A/ = l vector multiplet. 

Our aim is to study string solutions in a generic vacuum on the Higgs 
branch ()3.3|) . It is clear that this solution is formed by both massive electro- 
magnetic and scalar fields as well as by massless scalars living on the Higgs 
branch. Similar to the approach of ref. ^U] (reviewed above in Sect. 2) 
we would like to regularize the problem at hand in the infrared giving these 
massless scalars a small mass mi and then taking the limit mx — > 0. We will 
do it in the next subsection. 



3.2 Softly broken M=2 QED 

It turns out that we can lift the Higgs branch of M— 1 QED considered in the 
previous subsection giving the massless fields small masses without breaking 
M— 1 supersymmetry. To do so following ref. [TJ] let us consider A/ = 2 QED. 
This theory besides fields which enter M= 1 QED contains also neutral chiral 
field A which interacts with charged matter via superpotential 

W N=2 = ^=QAQ-^= V A, (3.7) 

where Q and Q are charge superfields and we also introduce the FI F-term 
proportional to the complex FI parameter 

77 = £i + i&- (3.8) 

Two FI F-term parameters £^2 together with one D-term parameter £3 form 
a triplet with respect to global SU(2) R symmetry of A/ = 2 QED. 

Let us break A/= 2 supersymmetry down to M— 1 adding a mass term for 
the chiral field A via additional superpotential 

SW = ^A 2 . (3.9) 

Now if we take the limit fi — > 00 we can integrate out heavy A-field and end 
up with M— 1 QED considered in the previous subsection. Namely, adding 
(13. 9j) to (j3.7j) and integrating out A we get the superpotential 

W = -^-(qQ-X) . (3.10) 



4/i V 2 
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This superpotential leads to the following scalar potential in A/= 1 QED 

mm 

V(q, q) = 9 j(\q\ 2 ~ \q? ~ &) a + ^{\q? + |sf " \? (3.11) 

The first term here comes from the D-components of the gauge multiplet, 
while the second one comes from the superpotential above. We expand rj = 
£1 + where £1 and £ 2 are real and use SU(2) R rotation to put £ 2 = 0. 

We see that in the limit \i —>■ oo this potential reduces to the one in 
(13. 2|) . However, for any finite fi the potential (j3.11|) has no Higgs branch 
and the vacuum state is uniquely defined. Namely, the potential ()3.11|) has 
a minimum at 

< q >= V« cosh Po, 



< q >= V&sinhpo, (3.12) 
where we introduce parameter p defined via 

sinh2p = 7- (3.13) 

Note, that if ^ 2 = we can use gauge freedom to make both < q > and 
< q > positive. 

Calculating the 4x4 scalar mass matrix near this vacuum we obtain one 
zero eigenvalue (corresponding to the state "eaten" by the Higgs mechanism) 
and another one equal to the mass of the photon, see ()3.4|) . ()3.6|) . As we 
explained above the corresponding scalar together with the photon form a 
bosonic part of 7V=1 vector massive multiplet. The mass of this multiplet 
remains unchanged by our IR regularization. The only modification is that 
the VEV of the scalar field is now fixed by the parameters of the theory. 
Namely, it is defined by ()3.5|) which can be now expressed in terms of FI 
parameters as 



v 



C 2 • C (3.14) 
The other two eigenvalues of the scalar mass matrix are given by 

m L = ^. (3.15) 

This is the mass of one A/=l chiral multiplet (containing two real scalars). 
In the limit of large fi that we consider here 

mi <^ m T (3.16) 
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This means that we have one heavy scalar in our problem ( with the mass 
equal to the mass of the photon, m# = m 7 ) and two light scalars. In 
particularly, in the limit of fi —>■ oo ttil goes to zero and we recover the Higgs 
branch of Af= 1 QED. 

Note that the form of the second term in the potential (j3.11J) is not 
important for our purposes. It serves as a IR regularization which lifts the 
Higgs branch and gives would be massless moduli fields a small mass (|3.15|) . 

4 String solution 

In this section we consider solutions for flux tubes which arise in a generic 
point on the Higgs branch of A/= 1 QED. First we start with a base point 
on the Higgs branch, which corresponds to < q >= 0. Next we consider the 
generic point on the vacuum manifold. 

4.1 BPS strings 

Consider the particular vacuum with < q >= in M— 1 QED (j3.1|) . It is well 
known that there is a BPS ANO string solution for this particular choice of 
vacuum [21]. 

We can easily recover this solution in the softly broken QED with the 
scalar potential (|3.11j) taking the limit £i = 0, while keeping v = £3 non-zero. 
It is clear that light fields do not play any role and we can look for the string 
solution using the following ansatz 

q = 0. (4.1) 

With this substitution the bosonic part of softly broken Af= 2 QED reduces 
to the standard Abelian Higgs model ()2.3|) for one complex scalar field q in- 
teracting with the electromagnetic field, at the particular value of the quartic 
coupling A = g 2 /8. This value of A ensures that the mass of the scalar field in 
(12. 3|) is equal to the mass of the photon, see ()3.6|) . As we already explained 
in the previous section this is a consequence of M— 1 supersymmetry. 

Thus our model is on the boundary separating superconductors of the 
I and II type. In this case vortices do not interact. It is well known that 
vanishing of the interaction at uih = vn^ can be explained by the BPS nature 
of the ANO strings. The ANO string satisfy the first order equations and 
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saturate the Bogomolny bound (221 • This bound follows from the following 
representation for the string tension T, 



T = 2^3 n + J d 2 x j ^F lk + | (|g| 2 - £ 



1 2 



3 J 



(4.2) 

Here indices i,k = 1,2 denote coordinates transverse to the axis of the vortex. 
The minimal value of the tension is reached when both positive terms in the 
integrand of Eq. ()4.2|) vanish. The string tension becomes 



T, 



BPS 



(4.3) 



where the winding number n counts the magnetic flux 2n n (we assume pos- 
itive n). The linear dependence of string tensions on n implies the absence 
of their interactions. 

For simplicity we consider the case n — 1. Putting the integrand in 
Eq. (J4.2j) to zero gives two first order differential equations, 



r— 0(r) - f{r)(p(r) = 
dr 

- £/(r) + £ (> 2 (r) - £ 







, , > ^ , , (4.4) 

r dr 4 v ' 

where the profile functions 0(r) and f(r) are introduced in a standard way, 



q[x) 
AAx) 



[r ) e 

x 



-2e l3 -± [1 - /(r)] 



(4.5) 



Here r = yxf is the distance and a is the polar angle in the (1,2) plane. The 
profile functions are real and satisfy the boundary conditions 



0(0) = , f(0) = 1 , 
0(oo) = v/& , /(oo) = 



(4.6) 



which ensures that the scalar field reaches its vev at the infinity and 
the vortex carries one unit of the magnetic flux. The equations (|4.4|) with 
boundary conditions ()4.6|) lead to the unique solution for the profile functions 
(although an analytic form of this solution is not found). 
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In QED with the M= 2 supersymmetry broken down to M= 1 the emer- 
gence of the first order equations (J4.4)) signals that some (half) of the remain- 
ing four SUSY charges of J\f— 1 algebra act trivially on the ANO solution (cf. 
[2lfl |2H |2TJ HI])- in this case the Bogomolny (topological) bound for the 
string tension coincides with the central charge of SUSY algebra. Note that 
at /i = we have BPS strings in Af=2 QED ITT] . As we increase /i but 
keep £i = they become BPS strings of M= 1 QED. As the number of states 
in string multiplets cannot jump we conclude that we get two BPS string 
multiplets at large \x from one BPS string in M=2 theory at // = 0. 

4.2 String in a generic vacuum 

Now let us turn to the generic case when both £i and £3 are non-zero. It is 
clear that in this case the ANO string is no longer BPS saturated. To see 
this, note that now we cannot take the light scalar fields with mass to 
be zero on the vortex solution. This means that the vortex has a long range 
"tail" in (1, 2) plane formed by the light scalars. This ensures attraction of 
different vortices and type I superconductivity. 

The string tension for non-BPS string is bounded from below by the 
central charge of the M= 1 algebra 

T > 2< 3 , (4.7) 

which is determined by the FI D-term parameter £3 |21j . 

Now let us work out the approximate solution for the ANO vortex and 
calculate its tension. The difference of the problem at hand with the one 
studied in [TU] and reviewed in Sect. 2 is that in [TU] was considered the case 
which corresponds to £3 = in the low energy QED. In this case heavy scalar 
does not develop VEV and can be put to zero on the vortex solution. The 
only scalar which is present in the problem is the light one. 

In the present case, when both ^ and £3 are non-zero both heavy and light 
scalars are non-zero. To solve the problem we adopt the following model for 
the structure of the vortex in (1, 2) plane. The electromagnetic field together 
with the heavy scalar form a core of relatively small radius. Let us call it 
R c . Outside this core heavy fields are almost zero, while light components 
produce a "tail' of size l/m^. To be more specific we consider separately 
three different regions: r < R c , R c <^r ^ l/m L and r > l/m L . 

In the first region (inside the core) 

r < R c (4.8) 
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only heavy scalar and electromagnetic fields are non-zero, while the light 
scalars are almost zero. This suggests that we can look for the vortex solution 
with 

q = 0. (4.9) 

With this ansatz our model reduces to the standard Abelian Higgs model 
()2.3|) at the special value of quartic coupling A = g 2 /8, see Sect 4.1. This 
means that the solution for the vortex inside the core is given by the BPS 
solution we reviewed in the previous subsection. In particular, the field q 
and Ai are given by ()4.5|) with the profile functions and / subject to the 
first order equations (|4.4|) and boundary conditions ()4.6|) . Now the bound- 
ary condition (J4.6j) at infinity should be understood as a condition at the 
boundary of the core , namely, at r > R c . Say, for scalar fields this means 



q(r > R c ) = ^3, (4.10) 

while q is zero. This means that outside the core scalar fields go to the base 
point on the vacuum manifold. 

As soon as the size of this BPS string is given by 1/ y/g 1 ^, we conclude 
that 

Rc = ^ r - (4.11) 

The string tension of the vortex with n = 1 is given by 

T = 2tt& + T tail . (4.12) 

Here the first term is "BPS" contribution ()4.3|) coming from the region ()4.8|) 
inside the core, while T ta n stands for the contribution of the light scalar "tail" 
coming from the region outside the core. Let us work out this contribution. 
Consider the region of intermediate r outside the core, 

R c <^r « l/ mi . (4.13) 

In this region we can neglect the second term in the potential (|3.11j) . The 
motion of scalar fields is restricted by the constraint 

\q\ 2 -\q\ 2 = & (4.14) 

enforced by the first term in the potential (|3.11|) . see also (|3.2j) . The con- 
straint (J4.14j) ensures that the only light scalar modes (actually massless, 
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once we neglect m L ) are exited in the region (|4.13|) . To put it another way, 
the constraint (|4.14j) reduces the motion of scalar fields to the motion along 
the Higgs branch given by ()4.14|) . From ()4.14|) we see that this Higgs branch 
has non-flat metric (see also eq. (j4.18|) below). This makes a difference with 
the case studied in ^Uj where Higgs branch was flat. 

Now let us work out the sigma model on the Higgs branch. The constraint 
(I4.14|) can be solved by the substitution 

q = ^~ 3 e ia+i(3(x) cosh p(r), 

~q=^ e fo -^W S inh p(r) . (4.15) 

Two complex fields q and q have two different phases. We use the gauge 
freedom to fix their average to be equal to the polar angle a in order to ensure 
the correct flux of the solution. The phase difference contains arbitrary 
function /3(x). The boundary condition for functions p and ft is fixed by 
VEV's of scalar fields (jXP2j) . 

p{r ~ R c ) = 0, 

p(r ~ l/m L ) = po, ft{r ~ l/m L ) = 0. (4.16) 

Now let us take the low energy limit formally sending m 7 — > oo. In this 
limit we can integrate out the gauge fied 

. qdjq - djqq + qdjq - dm 

A-i — —l ; — — 

qq + qq 

where we use the substitution (|4.15j) . Then our model reduces to the 2d 
sigma model 

Ttau = 6 / d 2 x cosh 2p {{d iP f + {diftf tanh 2 2p] (4.18) 

with non-flat metric of the target space. 

Now the problem is to find the classical solution for two dimensional sigma 
model (|4.18jl with boundary conditions (|4.1fi|L To do this in more general 
setting we consider a sigma model with arbitrary metric 

Ttau = £3/ <Pxg MN di(p N di(p N , (4.19) 
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where N, M numerates light scalar fields. Now we assume that these fields 
depend only on the radial coordinate r. This leads us to the one dimensional 
sigma model which determine the tension of the string "tail" 

T ta u = ; ^ = C dt g M Nd t <p N d t <p N , (4.20) 
log l/m L R c Jo 

where we introduce normalized logarithmic time 

log r I R c 

t=- — V 4.21 
log l/m L R c v ; 

The equations of motion for this sigma model define a geodesic line 

3 t V + T N MK d t <p M d t <p K = 0, (4.22) 

where denotes the connection on the target manifold. The energy con- 
servation shows that the action on this line is determined by the square of 
the length of this line 

= - 2 ^ 3 - I 2 . (4.23) 
logl/m Lj R c 1 ; 

Here the length of the geodesic line reads 



1 dt^g MN d t <p N d t <p N . (4.24) 



For our model (|4.18|) at hand the geodesic line is particularly simple. 
Clearly, 

(3 = (4.25) 
at the geodesic line and its length on the Higgs branch becomes 

I = £° dp^J cosh 2p, (4.26) 

where the upper limit is defined by f)3. 13|) . This length determine our final 
result for the tension of the string 

T = 2^ 3 + 2 ^_ 3 - I 2 . (4.27) 

log {gV&/m L ) 

It is easy to check that the third region r > l/rrii (where scalar fields ap- 
proach their VEV's exponentially) gives corrections to this result suppressed 
by powers of log g\/iz/m Ll see also [10 . 
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Let us note that this string solution "feels" not only the VEV of the scalar 
field v 2 but the whole structure of the Higgs branch. In particular, the size 
of the core R c is determined by g\f£$ (see (14.11(0 rather then by the mass of 
the photon (|3.4|) ( determined by gv). 

Of course, if we send £1 to zero going to the base point on the Higgs 
branch < q >= the second term in (J4.27j) vanishes and we will get the 
BPS string. Note however, that the long non-BPS string multiplet contains 
two bosonic and two fermionic states, while the short BPS multiplet contains 
one bosonic and one fermionic state. As the number of states cannot jump 
the long multiplet reduces to two short BPS multiplets at £1 = 0. This is in 
accordance with our conclusion in the end of the previous subsection where 
we considered breaking of M= 2 supersymmetry by turning on \x. There we 
saw that one A/=2 BPS multiplet reduces to two M= 1 BPS multiplets as we 
switch on jj, at zero £i. 

If £3 and v 2 are of the same order the second term in ()4.27|) becomes small 
as compared with the first one due to the logarithmic suppression. This is 
the reason why we can call this string "almost BPS". Note however, that 
this string is not a BPS one. It does not belong to a short BPS multiplet 
and all four supercharges act on this solution non-trivially. 

To make contact with the case discussed in Sect. 2 let us consider the 
limit ( 3 <C v 2 . In this limit the geometry of the Higgs branch becomes flat 
and I 2 — > v 2 /l; 3 . Then our result (|4.27|) reduces to eq. ([2.9)1 . Note that as 
we send £3 to zero the size of the core of the string R c grows and eventually 
freezes as it reaches the value R g , see ()2.7|) . At g\f£a <C ^/Rg the core is 
determined by the electromagnetic field JU] and the heavy scalar plays no 
role any longer. 

Let us now remove our IR regularization sending \i — > 00. Then we recover 
the Higgs branch of M= 1 QED and light scalars become strictly massless. 
Repeating the arguments in Sect. 2 we see that infinitely long strings do not 
exist any longer in the generic point on the Higgs branch. However, strings 
of finite length still exist. Their string tension is now given by 

T = 2< 3 + 1 J*% T J 2 - (4-28) 
log (0 Vf 3 L) 

They give rise to the confinement of monopoles with the potential of type 

OD. 

To conclude this section let us make a comment on literature. In ref. 
[ToT ITd] vortices in M= 1 QED were considered and it was concluded that in 
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generic point on the Higgs branch the string is unstable. The only vacuum 
which support string solutions is the base point of the Higgs branch < q >= 0. 
The so called " vacuum selection rule" was put forward in [HI] to ensure 
this property. We would like to stress that our results here and in ref. JU] do 
not contradict above mentioned papers. As we already explained infinitely 
long strings do not exist in a generic point on the modular space, indeed. 
However this does not mean that we loose confinement. As we explained in 
Sect. 2 to study confinement we have to consider strings of finite length. 
Strings of finite length do exist in a generic point on the Higgs branch and 
produce confining potential ()1.2|) . 

5 7V=2QCD 

In this section we study another example of flux tubes on Higgs branches. We 
consider A/ = 2 QCD with gauge group SU(2) and two flavors of fundamental 
matter (quarks). If masses of these quarks are equal then there is a Higgs 
branch in this theory. First we review the effective low energy theory on the 
Higgs branch [5 and then construct the string solution. 

5.1 Higgs branch 

The M= 2 vector multiplet of the theory at hand on the component level 
consists of the gauge field A^, two Weyl fermions Af a and A" (a = 1, 2) and 
the complex scalar <p a , where a = 1, 2, 3 is the color index. Fermions form a 
doublet A" a with respect to global SU(2)r group, / = 1,2. 

The scalar potential of this theory has a flat direction. The adjoint scalar 
field develop an arbitrary VEV along this direction breaking SU (2) gauge 
group down to U(l). We choose (ip a ) = 5 a3 (a). The complex parame- 
ter (a) parameterize the Coulomb branch. The low energy effective theory 
generically contains only the photon = A^ and its superpartners: two 
Weyl fermions Aj Q and the complex scalar a. This is massless short vector 
M—2 multiplet. It contains 4 boson + 4 fermion states. W-boson and its 
superpartners are massive with masses of order of (a). 

Quark hypermultiplets have the following structure. They consist of com- 
plex scalars q kA , q~Ak and fermions ip kaA , 4>A k , where k — 1, 2 is the color index 
and A = 1, . . . , Np is the flavor one. Scalars q kA , q kA form a doublet q k ^ A , 
f = 1, 2 with respect to SU(2)r group. All these states are in the BPS short 



17 



representations of M= 2 algebra on the Coulomb branch with 4x N c x Nf = 16 
real boson states (+ 16 fermion states). 

Coulomb branch has three singular points where monopoles , dyons or 
charges become massless. Two of them correspond to monopole and dyon 
singularities of the pure gauge theory. Their positions on the Coulomb branch 
are given by |3| 

u m4 = ± 2mA 2 - ^A 2 , (5.1) 

where u = \{^ a2 ) and A 2 is the scale of the theory with Nf = 2. In the 
large m limit u m ,d are approximately given by their values in the pure gauge 
theory u m ^ ~ ±2mA 2 = ±2 A 2 , where A is the scale of Nf = theory. 

The charge singularity corresponds to the point where half of quark states 
becomes massless. We denote them q^ A and ip aA , ip^ dropping the color 
index. They form Nf = 2 short hypermultiplets with 4 x Nf = 8 real boson 
states. The rest of quark states acquire large mass 2m and we ignore them 
in the low energy description. The charge singularity appears at the point 

a = - V2m (5.2) 

on the Coulomb branch. In terms of variable u ()5.2j) reads 

u c = m 2 + ~A\. (5.3) 

Strictly speaking, we have 2 + Nf = 4 singularities on the Coulomb branch. 
However, two of them coincides for the case of two flavors of matter with the 
same mass. 

The effective theory on the Coulomb branch near the charge singular- 
ity (J5.2)) is given by Af = 2 QED with light matter fields qf A , and their 
superpartners as well as the photon multiplet. 

We deform the underlying non-Abelian theory adding a superpotential 

5W = | $ a2 , (5.4) 

which is a mass term for the adjoint chiral field $ a . Generally speak- 
ing this perturbation breaks M= 2 supersymmetry down to M= 1 . The 
Coulomb branch shrinks to three above mentioned singular points which we 
call M= 1 vacua. In particularly, here we will be interested in quark vacuum 
(|5.3jl which is far away from the origin in week coupling provided the mass 
of quarks is large m ^> A 2 . 
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In the low energy effective theory near quark vacuum the superpotential 
(|5.4j) can be expanded as 



5W = -^= V 5A + ---, (5.5) 

where A is the neutral chiral superfield with the lowest component a, while 

rj = -2\/2 < a >= Afim. (5.6) 

Dots in ()5.5|) denote higher orders in 5 A = A + \/2m. 

It turns out that if in the limit of small /1 we truncate the series in ()5.5|) 
restricting ourselves only to the leading in SA term then the perturbation 
(15 .5j) does not break M= 2 supersymmetry in the low energy QED ■ 
The reason for this is that the leading term in (J5.5|) is linear in 5 A. It is FI 
F-term which does not break Af= 2 supersymmetry. 

The bosonic part of the effective Af=2 QED near quark vacuum looks 

like 

S QED = J ^xl-^F^ + l|0„ a | 2 + V^V^ 



g 2 



[ Tr (qr m q) - £ m ] 2 j + ^\q fA \ 2 \a + V2m A f 



(5.7) 



where trace is calculated over flavor and SU (2) R indices while £ m , m — 1,2,3 
is a SU (2)_r triplet of FI parameters. In particular, for the choice (j5.5J) ^3 = 
while £1 and £2 are real and imaginary parts of the complex parameter rj, see 
(J3.8j) . The scalar potential in the theory (J5.7)) comes from the elimination of 
D and F terms. In more transparent notations it reads 

2 2 c 

t/ 3 i\ A\2 |~ |2\2 1 9 I ~ A Si |2 

V = —{\q I - \qA\ ) + -jWAq ~ —\ 

+ \{\q A \ 2 + \q A \ 2 )\a + V2m A \ 2 . (5.8) 

The QED coupling constant g 2 is small near the quark vacuum in ()5.7j) if 
m > A 2 . 

The charge singularity ()5.2|) is the root of the Higgs branch 5j. To find 
it we look for zeros of the potential in (|5.7|) . We have 

qA P (r m ) p f q fA = f m , m = l,2,3. (5.9) 
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(Here m is an adjoint SU{2) R index, not to be confused with color indices.) 
This equation determines the Higgs branch (manifold with (q) ^ 0) which 
touches the Coulomb branch at the point ()5.2j) . It has non-trivial solutions 
for Nf > 2 5\. This is the reason why we choose Nf = 2 for our discussion. 

Once quark fields develop non-zero VEV's on the Higgs branch the U(l) 
gauge group in (J5.7)) is broken and the photon acquires the mass 



1 2 2 



where we introduce the quark VEV 

K<z M >l 2 H</> 



(5.10) 



(5.11) 



As soon as the potential is zero on the fields which satisfy constraint (|5.9|) 
the moduli fields which develop VEV's on the Higgs branch are massless. 

The number of these massless moduli (dimension of the Higgs branch) 
is four. To see this, note that we have eight real scalars subject to three 
conditions (|5.9jl . Also one phase is gauged. Overall we have 8-3-1=4, which 
gives us the dimension of the Higgs branch. Four massless scalars correspond 
to the lowest components of one short hypermultiplet. The other quark 
fields (4 real boson states + fermions) acquire the mass of the photon (J5.10)) . 
Together with states from the photon multiplet they form one long (non- 
BPS) Af=2 multiplet (cf. [H]). It has 8 boson + 8 fermion states. 

If we consider energies much less then the photon mass we can integrate 
out heavy scalars and electromagnetic field. Then we are left with the ef- 
fective sigma model for light fields living on the Higgs branch. The four 
dimensional Higgs branch is a hyper-Kahler manifold. At non-zero £ m it has 
Eguchi-Hanson geometry [T7J. Like in Sect. 3 we use SU(2)r rotations to 
put £2 = so 1] is real, rj = £1. The convenient parametrization of the metric 
is as follows pp] 



fd 4 x^ 


l-(% 2 




w z 



(d,w) 2 + w 2 



(9J) 



2 + sin 2 9{djf 



1 - 



6_x2 

2' 



W 



(5.12) 



Here the Higgs branch is parametrized by one modulus field w (which takes 
values from to 00) 

,M|2 



w 



Vf 



\<1 A \ 2 + \U\ 2 



(5.13) 
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and three phases 9 ( with values in the interval (0, 7r)) and 5, 7 with values 
in the interval (0, 2ir). 

We use three rotations of broken ST/ (2) subgroup of global SU(Nf = 
2) x SU(2)r group to put VEV's of the scalar fields on the Higgs branch in 
the form 

<w>=Vi < e >= 0, < 7 >= 0. (5.14) 

One important distinction of the Higgs branch at hand with the one in 
J\f= 1 QED is that the base of the Higgs branch in the present case is a 
compact manifold rather then a point. The base of the Higgs branch is 
defined as a submanifold with the minimal \q\ 2 . For the case of Higgs branch 
(13. 3j) in M— 1 QED the base is defined by the condition q = which reduces 
the Higgs branch to a single point p = 0. For the case of the Higgs branch 
(15. 9j) this condition becomes 

q A = q A , (5.15) 

which can be obtained from the condition q = by SU(2)r rotation trans- 
forming £3 into £1, see [11J. The condition ()5.15|) reduces the number of real 
scalars from eight to four. They are subject to constraint ()5.9|) which boils 
down to a single condition 2|g A | 2 = ^. Subtracting U{1) phase we get 4- 
1-1=2 which is the dimension of the base of the Higgs branch. Clearly this 
manifold is a two dimensional sphere S2. In terms of the coordinates in (|5.12|) 
it is parametrized by angles 9 and 5, while 

w = ^, (5.16) 

on the base manifold. 
5.2 Flux tubes 

Now let us consider ANO strings in a generic point on the Higgs branch in 
W=2 QED ipT7|) . The field a is frozen at its VEV 

a = —\pm (5-17) 

and does not play any role in the string solution. As we discussed in the 
previous subsection we have massive scalars with the mass equal to the mass 
of the photon (|5.10j) and four massless scalars. Therefore to find the string 
solution we use the same method as in Sect. 4. 
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Namely, our string consist of a BPS core formed by heavy fields and a 
tail formed by light fields. To find the BPS core we impose condition 



q A = r = ^ A , (5.18) 

which reduces the QED (15. 7|) to the Abelian Higgs model with two complex 
flavors ip A with the potential 

9 2 (\A\1 



Clearly this model possess standard BPS strings 1 with the tension 

T BPS = 27T& (5.20) 
which satisfy boundary conditions 

w(r > R c ) = yfe, 6(r > R c ) = 6 , 5(r > R c ) = S , (5.21) 
outside the core. Here the size of the core R c is given by 

Rc = -^=. (5.22) 

Outside the core heavy fields are almost zero and the string is determined 
by the classical solution of the one dimensional sigma model (|4.2U|) (with £ 3 
replaced by £i) with target space geometry ()5.12j) . Namely, the tension of 
the tail is given by 

2ir 2 

Ttaii = \ pr^- I , (5.23) 

logL/i? c 

where L is the length of the string and I is the length of the geodesic line on 
the Higgs branch (|5.12|) between points (|5.21|) and (|5.14|) . 

The initial point of this geodesic line ()5.21|) is not fixed yet. In principle, 
it could be any point on the base submanifold. To fix it we require the string 
to have the minimal tension. This boils down to the condition of having 
geodesic line of the minimal length between point (J5.14j) and the point in 
question on the base submanifold. 



^ee however next section for the discussion on semilocal strings 
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This requirement ensures that this point has 



w(r > R c ) = ^l, 9{r > R c ) = 0. (5.24) 

Moreover, clearly the phases 9, 5 and 7 are zero on the whole geodesic tra- 
jectory. We checked this explicitly solving the Hamilton- Jacobi equation for 
this geodesic line following the method of ref. |26j . This means that the 
length of the geodesic line is given by 

I = / ^ , 5.25 

Ji-(h) 2 



With this length the final answer for the tension of the string becomes 

T = 2^ 1+ I 2 . (5.26) 

log (#V6i L ) 

If we take the VEV's of scalar fields on the base submanifold (v = v^T) 
the second term in (|5.26|) vanishes and we get tension of BPS string. Note, 
that like in Sect. 4 we get two short BPS multiplets in this limit from one 
long non-BPS string multiplet which exist in a generic point on the Higgs 
branch. In the opposite limit £1 —>■ the Eguchi-Hanson manifold becomes 
flat and / — > v. Besides that the BPS core contribution in (J5.26j) vanishes 
and we get (|2.12|) obtained in [TU] for the case of Higgs branch without FI 
term 2 . 

To conclude this section we can reexpress the boundary conditions on the 
base manifold ()5.24j) which serves as a conditions "at infinity" for the BPS 
core of the string in terms of the original quark fields (p A of the model with 
potential (|5.19|) . Using the standard relation of unit vector on S 2 in 0(3) 
sigma model with quark variables 

cos# = -Cpa{t 3 )bV B , 
Si 

sin 9 cos 5 = -^<Pa(t 1 )b 1 P B \ 
4i 

sin 9 sin 5 = ^ a (t 2 )^ b (5.27) 
£1 



2 Note that the size of electromagnetic core in the limit £1 — > is frozen on the value 
R g , see 
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we obtain that point (|5.24j) on the base submanifold corresponds to the fol- 
lowing values of quark fields 

V A=1 = s/li, V A=2 = 0. (5.28) 

This means that with our choice of scalar VEV's ()5.14|) the BPS core is 
formed only by the first flavor, while the second one remains unexcited. 

6 Semilocal strings 

So far our study of flux tubes in M—2 QED was not complete. The point 
is that BPS strings in theories with extended global symmetry possess an 
additional zero mode associated with their transverse size r$. The string 
parametrized by the size r is called semilocal string (see ^H] for a review). 
Semilocal string interpolates between ANO string and 2d sigma model in- 
stanton lifted in four dimensions (lump). At non-zero r the semilocal string 
has power fall-off of the profile functions at infinity, instead of the exponential 
fall-off for ANO string at r = 0. 

As we explain below this leads to a dramatic physical effect - we loose 
confinement if a semilocal string is developed instead of ANO string. In 
the next subsection we briefly review semilocal strings and then consider 
the possibility of their formation in a generic point on the Higgs branch in 
M— 2 QCD with non-zero FI term. 

6.1 BPS semilocal strings 

Let us recall basic features of semilocal strings [18j . The simplest model 
where they appear is the Abelian Higgs model with two complex flavors 

Sah = /^|^Fi+|V^| 2 + ^(|^| 2 -6) 2 }. (6.1) 

Note, that we already considered this model in Sect. 5.2 discussing the BPS 
core of the string on the Higgs branch. It arises from M— 2 QED ()5.7j) when 
we restrict ourselves to the ansatz (|5.18j) . 

The topological reason for existence of ANO vortices is that for gauge 
group U(l) 7Ti[U(l)] = Z. On the other hand we can go to the low en- 
ergy limit in (j6.1|) restricting ourselves to the vacuum manifod \tp A \ 2 = £i. 
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The vacuum manifold has dimension 4-1-1=2, where we subtract one real 
condition mentioned above as well as one gauge phase. It represent two di- 
mensional sphere 5*2. Thus, the low energy limit of theory ()6.1|) is 0(3) sigma 
model. Now recall that ^[S^] = 7i"i[T/(l)] = Z and this is a topological rea- 
son for existence of instantons in two dimensional 0(3) sigma model. Lifted 
in four dimensions they become a string-like objects (lumps). 

So now the question is what is the relation between ANO flux tubes of 
QED ()6.1|) and lumps of 0(3) sigma model. Clearly the model ()6.1|) has 
ANO strings. Say, if we put the second flavor to zero this model reduces to 
the one flavor model considered in Sect. 4.1. Then ANO string is given by 
eqs. (|4.5j) . (|4.4|) for the first flavor while the second one is zero. However, it 
turns out (see [27]) that this solution has zero mode associated with exiting 
of the second flavor. This zero mode is parametrized by the parameter ro 
which plays the role of the size of the string in (l,2)-plane. 

To find this solution let us modify the standard parametrization ()4.5|) for 
the ANO string including the second flavor 



ip 2 (x) = fa(r), 

A(x) = -2ey ^ [1 - f(r)) . (6.2) 

Note, that the second flavor does not wind at infinity. Therefore, it has 
boundary condition 2 (oo) = while at r = it can be non-zero. 
The first order equations for the profile functions here look like 

r^-0i(r) - /(r)0i(r) = , 

r^Mr)-(f(r)-l)(h(r) = , 

- \ ^f(r) + 9 j (#(r) + $(r) - £i) = . (6.3) 

The solution to these equations fI7\ |2*H] at non-zero r is very different from 
that of ANO string. It has long range power fall-off at infinity for all profile 
functions. In particularly, in the limit of large transverse size of the string 
r l/<?\/Ci ^ nas the form 



Mr) = V£ 



r 2 + Tn 
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02 (r) 



2 _i_ r.2 ' 
"I" ' n 



f = -T^- ( 6 - 4 ) 
This solution has the same tension as ANO string 

T = 2vr£i. (6.5) 

We see that scalar fields on the solution ()6.4|) belong to the vacuum 
manifold \f A \ 2 = £i. This means that we can relate this solution to the 
0(3) sigma model lump. To do this we first use the relation between fields 
tp A and the unit vector on S2 f|5.27|) and then represent this vector in terms 
of complex field u(x) via standard relations 



cos 9 



2 



sin 9 cos 5 = 2 



1 + M 2 ' 

Reu 



1 + M 2 ' 



Itti uj 

sin 9 sin 5 = 2 : — --. (6.6) 

1 + \uj\ 2 

With this substitution the low energy limit of (j6.1j) becomes an 0(3) sigma 
model 

S W A (TW (6 - 7) 

The standard lump (instanton) solution of this model with center at zero 
and size r$ looks like 

Ulump = — , (6-8) 

Xi + IX 2 

where parameter r can be made real by the shift in the polar angle a. 

If we now reexpress this solution in terms of quark fields using relations 
(|6.6p . (|5.27J) we arrive at the solution ()6.4j) . Thus we identified the semilocal 
string in the limit of large tq (16. 4|) with the lump solution of 0(3) sigma 
model (up to a factorization over Z2). The reason for this identification is 
that 0(3) sigma model is a low energy effective theory for two-flavor Abelian 
Higgs model (j6.1|) . 
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Now let us come back to our Af=2 QED (|5.7|) and consider first the 
vacuum which belongs to the base submanifold v = a/6" of the Higgs branch. 
This was done in ^SJ^]. In this case the ANO string considered in the Sect. 
5.2 becomes BPS. As we already mentioned for these VEV's we can look for 
string solution using the ansatz (I5.18J1 which reduces the bosonic part of our 
theory to the two flavor model (jbMj) . This model has semilocal strings so 
our BPS ANO string is just a particular case of a semilocal string at r = 
[T3J[Tn|. Say, in the opposite limit r 3> the semilocal string is given 

by the profile functions (J6.4|) . 

The most physically important consequence of emergence of semilocal 
strings is that we loose the monopole confinement. As we already explained 
to study confinement we have to consider a string of a finite length L stretched 
between heavy monopole and anti-monopole. Clearly the problem now be- 
comes three dimensional. The string size r does not correspond to a zero 
mode any longer. Instead, it is fixed by the separation L at large value 
tq ~ L. Thus the monopole flux is not trapped into a narrow flux tube. In- 
stead, it is spread over a large three dimensional volume of size of order of L. 
Clearly this produces a potential between monopole and anti-monopole with 
power-like fall-off at large separations, V(L) ~ L~ 7 (7 > 0). Such potential 
does not confine. 

We see that formation of semilocal strings at a base point on the Higgs 
branch lead to a dramatic physical effect - deconfinement. Therefore it is 
quite important to study the possibility of formation of semilocal strings at 
a generic point on the Higgs branch. We do it in the next subsection. 

6.2 Semilocal strings at a generic vacuum on the Higgs 
branch 

Now we assume that semilocal string can be formed at a generic vacuum on 
the Higgs branch of Af=2 QCD and study the dependence of its tension on 
its size tq. As we mentioned in the previous subsection there is no such de- 
pendence for BPS string in the vacuum which belongs to a base submanifold, 
ro is associated with the zero mode. We will see below that this is not the 
case for a semilocal string at a generic vacuum on the Higgs branch. 

To study the possibility of formation of semilocal strings it is sufficient 
to consider a semilocal string in the limit of large size r 3> I/^a/^T because 
in this limit the semilocal string shows its crucial distinctions from the ANO 
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one. As we discussed in the previous subsection in this limit the semilocal 
string is formed by light fields only. It does not have core formed by heavy 
fields and looks like sigma model lump. 

Therefore to study semilocal string in the limit 3> 1/g v^T we can go 
to the low energy limit in J\f—2 QED ()5.7|) which is given by sigma model 
(I5.12j) with Eguchi-Hanson geometry. Assuming that all fields in this sigma 
model depend only on coordinates in (1, 2)-plane we rewrite the tension of 
our semilocal string as follows 



T. 



emilocal 



S d2x \ 


1 - (%' 




w z 



-I -1 



id,,™ 



) 2 + w 2 



{d lx 6f + sin 2 6{d^Y 



6_\2 
2> 



' W 



(<9 m7 + cos Od^f 



(6.9) 



Now our aim is to find a classical solution of this model which represents 
a kind of hybrid of the "tail" solution of Sect. 5.2 with lump (instanton) 
solution of Sect. 6.1. The boundary conditions of this solution at infinity are 
given by the VEV's (IH~TD 



w(oo) = v, 9(oo) = 0, 7(00) = 0, 



(6.10) 



while the boundary conditions at zero are given by those of the lump (J6.4J) . 
namely ^(0) = and y? 2 (0) = Rewriting this in terms of fields entering 
the model ()6.9j) using (|5.27|) we get 

(6.11) 



w(0) = ^ 6(0) = 7T, 7 (0) = 0. 

We see that our solution winds around base cycle of the Eguchi-Hanson 
manifold (S2) and then goes away in a non-compact direction (along w). 
Awaiting for the full solution of the problem here in this section we consider 
the case when the VEV of scalar fields is much larger then the FI parameter, 

v > V^- ( 6 - 12 ) 

In this limit the field w goes far away from the base (at w = y^T) along 
the non-compact manifold. When w 3> y^i the geometry in ()6.9|) becomes 
flat and the problem reduces to solving of free equations of motion. For the 
radial coordinate w the solution of the equation d 2 w = reads (cf. (|2.8|l ) 

log r/r Q 



w(r) = v 



log L/Tq 



(6.13) 
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while 

9{x) = 0, 8(x) = 0, 7(ar) = (6.14) 

like for the "tail" solution of Sect. 5. 2. Here we use the lump size r as a UV 
cutoff for the logarithmic behavior of w. This solution is valid at large r, 
R <C r <C L, where parameter R will be determined shortly. 

Instead for small r, r < R the solution is given by a certain deformation 
of the instanton solution of Sect. 6.1 which we denote as 

(w(x),6(x),6(x),j(x)) = {w inst (x),9 inst (x),6 inst (x),^ inst (x)), (6.15) 

where Wi nst (x) ~ Clearly the solution changers its behavior from ()6.15j) 
to ([6. 13)1 . (|6.14|) when the coordinate w given by the logarithmic expression 
f!6.13|) becomes much larger then y^T- This gives 

R = r exp (— ^ log Lira). (6.16) 

v 

The tension of our semilocal string is now given by the sum of tensions 
of deformed instanton ()6.15|) and the logarithmic "tail" ()6.13j) . ()6.14|) . 

2ttv 2 

T semilocal = Tinst[ r o) + "j T~l ' (6-17) 

logL/ro 

where the calculation of the tail contribution is similar to that in Sect. 2 (see 

Of course we do not know the function Ti nst {r§) without knowing the 
solution ()6.15j) . What we know is that at r = the semilocal string reduces 
to ANO string and deformed instanton ()6.15j) goes into BPS core of the string 
discussed in Sect. 5.2. Thus, 

T inst {r = 0) = 27r£i, (6.18) 

see (I5.20JI . Moreover, the central charge of M=2 algebra gives us a lower 
bound for this tension (see, for example, jTTj for details) 3 

T mst {r ) > 2<i. (6.19) 



3 To prove this formally we can construct configuration given by l|6.15|l at r < R and 
extrapolated by constant fields at r ^> R. Then we use the SUSY bound for this configu- 
ration. 
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We see that at least the second term in (jbM7|) depend on r , so r is not 
associated with conformal zero mode any longer. The reason for this breaking 
of conformal invariance at the classical level is the non-compactness of the 
Higgs branch. To see this note, that in the last subsection we were dealing 
with lump solution which maps two dimensional space onto compact base 
submanifold S2 of the Eguchi-Hanson Higgs branch. The tension of this lump 
solution is a constant (J6.5|) determined by the volume of S2. The collective 
coordinate ro in this case is associated with the conformal zero mode. 

As we move to a generic vacuum on the Higgs branch our semilocal string 
solution becomes a map onto a non-compact manifold. In order to get a finite 
string tension we use a cutoff in the target space introduced by the VEV of 
the scalar fields v. This cutoff in the target space requires an IR cutoff L 
in the two dimensional (1, 2)-plane. We already discussed how this works 
in Sect. 2. The logarithmic solution to the free equations of motion (j6.13|) 
cannot go to its VEV at infinity. In order to insure that the scalar fields 
reach their VEV's at infinity we have to introduce a IR cutoff, say small 
masses for light fields (like in Sect. 2,4) or the finite length of the string L. 
Clearly this breaks the conformal invariance and produces the dependence 
of the string tension (|6.17|) on ro- 

Now minimizing the string tension (|6.17|) with respect to ro using (|6.18|) 
and ()6.19|) we find that 

r = 0. (6.20) 

In the limit r — > the semilocal string reduces to ANO string considered 
in Sect. 5.2. In particular, the string tension is given by ()5.26|) in the limit 
I — > v. Note, that quantum fluctuations of r are negligibly small, suppressed 
by the world sheet volume of the string. 

This means that semilocal strings are not formed in the generic vacuum 
on the Higgs branch at least in the limit v ^> This conclusion is in 

accordance with the general expectation JH] that semilocal strings are not 
stable in type I superconductor and reduce to ANO strings. 

6.3 Deconfinement phase transition 

The conclusion made in the previous subsection that semilocal strings are not 
developed at a generic vacuum on the Higgs branch (at least if v ^> y^Y ) and 
we are dealing with ANO strings has rather important physical consequences. 
This means that we have a confinement phase for monopoles in these vacua 
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with confining potential determined by the ANO string tension (J5.26|) . see 
(|1.2|) . Instead, for the vacua on the base submanifold v = >/£i semilocal 
strings are formed and we have deconfinement phase. 

Therefore, it is clear that we can expect confinement /deconfinement phase 
transition at certain intermediate critical value v c . To find this critical value 
we have to study semilocal string solution for vacua which are not far away 
from the base submanifold, v ~ This is left for a future work. Here we 
put forward a plausible conjecture that in fact 



which means that we have confinement phase for all vacua on the Higgs 
branch except the base submanifold. The reason for this conjecture is that 
nothing special happens at any intermediate value v c . Instead, the value 
v c = is physically distinct because at this value of scalar VEV strings 
become BPS-saturated (one long non-BPS string multiplet becomes two short 
BPS multiplets, see Sect. 5.2). 

This perfectly matches results of jTU] where unbroken Af=2 QCD with 
two degenerative flavors at £1 = was considered. At zero £1 the geometry 
of Higgs branch becomes flat, given by R 4 /Z 2 and the S 2 cycle shrinks to 
a singular point at the origin. It was shown in ^U] that monopoles are in 
the confinement phase on the Higgs branch at any non-zero v. As v goes to 
zero the transverse size of ANO string becomes infinite while its tension goes 
to zero, see (j2.11j) . ()2.12j) . so monopole confinement becomes unobservable. 
Introducing the non-zero FI parameter £1 we resolve the singularity at the 
origin and see that at the base cycle S 2 we have deconfinement phase while 
on the rest of the Higgs branch we have monopole confinement. 

To conclude this section we would like to stress that the above discus- 
sion of confinement /deconfinement phase transition refers only to A/ = 2 QED 
(|5.7|) considered on its own right. In contrast, if we start from non-Abelian 
softly broken A/ = 2 SU(2) QCD the story is rather different. The point is 
that ( as we mentioned in Sect. 5.1) the Af=2 QED ()5.7|) is the low energy 
description of the underlying non-Abelian theory only in the special limit 



which ensures that M= 2 supersymmetry is preserved. In this case the su- 
perconductivity on the base submanifold is on the border between type I and 




(6.21) 



ji — > 0, m — > 00, £1 = const, 



(6.22) 
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type II and we have BPS strings which eventually become semilocal strings 
in the presence of SU (2) global flavor symmetry. 

Generically if we take into account corrections m/j/m the M= 2 supersymmetry 
is broken down to M= 1 . As it is shown in ^T] superconductivity becomes 
of type I. For the type I superconductor semilocal strings become unstable 
|18j . In particular, we expect that the string tension even for the string in 
the vacuum on the base submanifold becomes a function of r 

Tiump ~ 27r£i (l + const rl(5m) 2 H ), (6.23) 

where 5m is the splitting of masses in the former M=2 multiplet. This 
splitting is given by [TT] 

5m ~ u. (6.24) 

We conclude that in this case the semilocal string is not developed. Mini- 
mization of (|6.23|) with respect to r gives 

r = 0. (6.25) 

Note again that quantum fluctuation in r are suppressed by the volume of 
the string world sheet. 

This means that we do not have deconfinement phase even on the base 
submanifold of the Higgs branch in QCD with 7V = 2 supersymmetry softly 
broken down to M— 1 . All vacua on the Higgs branch of this theory give rise 
to the confinement phase for monopoles. 

7 Conclusions 

In this paper we studied flux tubes on Higgs branches in supersymmetric 
gauge theories. Our main conclusion is that although Higgs branches ensure 
presence of massless scalars flux tubes still exist and give rise to a confining 
potential. However due to the presence of massless scalars the string tension 
is no longer a constant. It becomes a slow (logarithmic) function of the length 
of the string L. Still it produces a confining potential of type (jl.2|) . 

First we studied M= 1 QED with FI .D-term. This theory has a two 
dimensional Higgs branch. The solution for the string in a generic point on 
this Higgs branch is given by "BPS core" formed by heavy fields and a "tail" 
formed by light fields living on the Higgs branch. Finding the solution for 
the "tail" reduces to finding a classical solutions of the sigma model on the 
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Higgs branch. The tension of the string is given by the sum of tensions of 
"BPS core" and the "tail", see §£231 . 

Next we studied M= 2 QCD with gauge group SU(2) and two flavors of 
quarks with the same mass. This theory has a four dimensional Higgs branch. 
We deformed this theory with the mass term for the adjoint field. To the 
leading order in the deformation parameter fi M— 2 supersymmetry is not 
broken in the effective QED describing the low energy limit of this theory. 
Higgs branch has an Eguchi-Hanson geometry. 

We showed that far away from the base cycle S2 of the Higgs branch 
stable ANO strings are formed. Their tension is again given by the sum 
of tensions of "BPS core" and the "tail", see ([5.26)1 . These strings produce 
confinement of monopoles with the confining potential of type ()1.2jl . How- 
ever, in the vacua on the base cycle 52 of the Higgs branch the story is 
rather different. Semilocal BPS strings are developed in vacua on the base 
submanifold leading to a deconfinement regime. We conjecture that the con- 
finement /deconfinement phase transition occurs exactly at the base cycle of 
the Eguchi-Hanson space. 

If we introduce breaking of M= 2 supersymmetry down to M— 1 in this 
theory it turns out that the deconfinement phase disappears and we have 
monopole confinement on the whole Higgs branch. 
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